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Abstract. We present a discrete multiscale expansion of the lattice potential Korteweg-de 
Vries (lpKdV) equation on functions of infinite order of slow-varyness. To do so we introduce 
a formal expansion of the shift operator on many lattices holding at all orders. The lowest 
secularity condition from the expansion of the lpKdV equation gives a nonlinear lattice equation, 
depending on shifts of all orders, of the form of the nonlinear Schrodinger (NLS) equation. 



1. Introduction 

Reductive perturbation techniques [14, 15] have proved to be important tools to find approx- 
imate solutions for many physical problems by reducing a given nonlinear partial differential 
equation to a simpler equation, which is often integrable [4]. Recently a few attempts to carry 
over this approach to partial difference equations have been proposed [3,6-9]. 

The basic tool of the discrete reductive perturbation technique developed in [7, 9] is a proper 
multiscale expansion carried out by introducing various scales and lattices defined on them. Let 
us recall some of the main results. The transformation between two different lattices of indices n 
and ni is given by [5] 

A>„ = g^rP v^^lA;^, (1.1) 

i=0 ^ ' i=j 

where u n : Z — ► R is a function defined on a lattice of index n G Z and u ni : Z — » R is the same 
function on a lattice of index m € Z. By the symbol A we mean the standard forward difference 
of the function u w.r.t. its subscript, e.g. A n u n = u n+ \ — u n . Here the coefficients P, j are given 
by 

k=j 

where u is the ratio of the increment in the lattice of variable n\ with respect to that of variable n. 
The coefficients Sf and & k are the Stirling numbers of the first and second kind respectively [2]. 
Eq. (jl.ip implies that a finite difference in the discrete variable n depends on an infinite number 
of differences on the variable n±, e.g. the function ii n +i can be written as a combination of the 
functions Uj's for i varying on an infinite subset of the lattice with index ni. Formula (II. 1|) can 
be inverted yielding 

i 

k=j 

In [7,9] one has introduced the notion of slow-varyness of order I for a function u n iff A^ +1 u n = 
(or equivalently A^ 1 u ni = 0, see [9] for further details). This definition enables to reduce 
infinite series to a finite number of terms. Moreover one has considered a generalization of 
formula (jl.ip in order to deal with functions u n = u n .^ n .y< depending on a finite number K 

of lattice variables rtj, 1 < i < K. The computations done in [7] proved that an integrable 

l 



A3 u - j\ V ^A' u 

i=j 
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lattice equation, as the lattice potential KdV (lpKdV) equation, reduces to a completely discrete 
nonlinear Schrodinger (dNLS) equation of the form 

i(0n,m+l - <fin,m) + Cl(0 n +l, m ~ 20 n ,m + <j>n-l,m) + C2<f>n,m\(t>n,m\ 2 = 0, (1.2) 

being c±,C2 two real coefficients. Let us stress the fact that the above dNLS equation has a 
completely local nonlinear part. We refer to the papers [7, 9] for further details. It has been 
proved by singularity confinement [13] and algebraic entropy [16] that the constructed dNLS 
equation is not integrable. 

In [4] Calogero and Eckhaus have shown that a necessary condition for the integrability of a 
nonlinear partial differential equation is that its multiscale reduction be integrable. Here, trying 
to find a lattice analogue of the Calogero-Eckhaus theorem, we extend the techniques developed 
in [7, 9] to any order of slow-varyness. 

In Section [2] we introduce a formal multiscale expansion holding at all orders of slow-varyness. 
Then in Section [3] we apply this technique to the lpKdV equation, thus providing an extension 
to all orders of the dNLS equation obtained in [7,9]. Finally, Section 0] is devoted to concluding 
remarks. 

2. Multiscale expansion on a lattice 

2.1. Lattices and shifts defined on them. Let u n : Z — > R be a function defined on a lattice 
of index n € Z. We can always extend it to a real function u : R — > R by defining a real 
continuous variable x = na x , where a x £ R is the lattice spacing. 

Let us define the shift operator T n such that T n u n = u n+ \. For u(x) we can introduce 
the operator T x , corresponding to T n , such that T x u(x) = u(x + a x ). The Taylor expansion 
of u(x + cr x ) centered in x reads 

2 i 00 i 

T x u{x) = u{x) + a x u^\x) + ^ 2 \x) + ... + ^-u^{x) + ... = V ?$u®(x), (2.1) 

i=0 

where 1*^(2;) = d l u(x)/dx l = d x u(x), being d x the total derivative operator. Eq. (12. ip suggests 
the following formal expansion for the differential operator T x : 



00 j 

i=0 



Introducing a formal derivative with respect to the index n, say S n , we can define, by analogy 
with T x , the operator T n as 



T n = e 5 -=J2^r- ( 2 - 2 ) 
i=o l - 

The formal expansion (j2.2[) can be inverted, yielding 

5 n = lnT n = ln(l + A n ) = V A^, (2.3) 

i=i 

where A n = A+ = T n — 1 is the discrete first right difference operator w.r.t. the variable n, see 
eq. (jl.ip . Notice that this is just one of the possible inversion formulas for the operator 5 n . For 
example it can be also written in terms of left difference operators A~ = 1 — T^ 1 : 



5 n = -lnT- 1 = -ln(l-A-) = J2 



(A- 



1 

i=i 



or in terms of symmetric difference operators A* = (T n — T n 1 )/2: 

5 n = sinh- 1 A^, = f;^M(A^, 



i=l 
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where Pi{x) is the i-th Legendre polynomial evaluated in x = 0. Hence the 5 n operators are 
formal series containing infinite powers of A n , but, acting on slow-varying functions of order £, 
they reduce to polynomials in A n of order at most I. Consequently, any formula written in terms 
of powers of S n for a given slow-varyness order £, contains its version for smaller orders j < £. 

The convergence of the series 5 n u n depends on the analyticity properties of the function u n . 
Hence, from now on, we will proceed formally, considering that the developments can be justified 
a posteriori, because the functions which finally appear (solutions to the resulting difference 
equations) will have the correct analyticity properties. 

2.2. Dilations on the lattice. Let us introduce a second lattice, obtained from the first one 
by a dilation. At first it is convenient to visualize the problem as a change of variable between 
the continuous variable i£K and a new continuous variable x\ = ex, < e <C 1. On the lattice 
one considers a change from the index n = x/a x to the new index n\ = xi/a Xl , where a xi is 
the new spacing. Assuming that a Xl ^> a x we can set a x = ea Xl , < e <C 1, so that n\ = een. 
As n,ni £ Z, ee is a rational number and one can define in all generality ee = M\/N <C 1 
with M±,N E N. However, if we want that the lattice of index n\ is a sublattice of the lattice 
of index n, we have also to require that M\/N = 1/M with M G N. 

The relationship between the discrete derivatives defined in the two lattices is given by 
eq. (jl.ip , for which we shall give a straightforward proof based on the well-known formulas [2] 



(e* - iy = j\ -^x\ [ln(l + x)] k =k\J2 -f x '' 

k=j i=k 

Here n\ = n{M\/N) and thus uj = M\/N. From the expression of An, j € N, we get by a 
straightforward algebra: 

00 ( A \k 

k=j 

k 00 t^j / \ k 



k=j k=j 



A 1 



k=j i=k i=] 

As a x = ea Xl and ee = M\/N eq. is proven. 

2.3. Discrete multiscale expansion. We present here the formulas necessary to construct 
a discrete multiscale expansion. According to the definitions given above, let us consider u = 
Un-m = u(x; x\) as a function depending on a fast index n and a slow index m = n(M\/N). At 
the continuous level, the total derivative d x acting on functions u(x;xi) is the sum of partial 
derivatives, i.e. d x = d x + ed Xl . What is the situation on the lattice? Let us construct a 
relation between the total shift operator T n and the partial shift operators T n ,T ni (T n u n - ni = 

H j n J c\\n\iTn\U"n\n\ — ^n;ni+l)- As 

— g^xdx _ g&x@Xg£&xOxi 



we can write 
with 



Tn = e Sn e (M 1 /N)8 ni ^ T n T^' h/N \ (2.4) 
T n = ± 6 f, Ti*M = £ , (2.5) 

i=0 ' i=0 

where 5 ni is given by eq. (|2.3|) with n substituted by m. 
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Eq. (|2.4p can be easily extended to the case of K slow variables = e l x, 1 < i < K . The 
action of the shift operator T n on a function u = u t \k depending on both fast and slow 
variables can be written in terms of the partial shifts T n ,T ni : 



K 



T n = T n \{T^\ (2.6) 



where the e ni 's are suitable functions of e and e depending parametrically on some integer 
coefficients Mj G N, 1 < i < K. 

To develop the fields appearing in partial difference equations with two independent discrete 
variables one has to consider the action of the operator (12. 6p on a function depending on two 
fast indices n and m, and on a set of K n + K m slow variables {ni}^ and {rrii} i J^, i.e. on 
u = u r x K n t \K m . Notice that in principle it is possible to consider K n = K m = oo. 
For the moment we assume a common definition of the small parameter e for both discrete 
variables n and m, but we denote with Mj the integers for the slow variables n« and with Mj 
the ones for mj. We have: 

Mi Mi 

In the next Section we shall consider a partial difference equation defined on a quadrilateral 
lattice, namely an equation of the type f(u, T n u, T m u, T n T m u) = 0. Assuming K n = 1 and K m = 
K we get from eqs. (12.51 12.61) the following expansions for the shift operators appearing in the 
equation f(u,T n u,T m u,T n T m u) = 0: 



1 n 1 n\ — 1 n 



(2.7) 



K 



t TTt 

J. J. 



(emj 



1=1 



T 

J. rn 



I 1 ( M 2 \ 

1 + ^M 1( 5 mi + ^ ( + M 2 6 m2 \ + 



N 2 2 



1 / M 3 \ 

+ W ~6~ 6 '^ + M i M 2^^ m2 + M 3 5 m:i + 0(1/N 4 



rp rp _ q-q-( e ni)q- I I q~( e 

- L n- L m — I n I ni I m\\ I mi 



K 



-( £ m,) 



T T 



1 + ^ (Mi<5 ni + Mi<5 mi ) + 



1 (M 2 2 



M 



+ N 2 ^ 2 ^ + M ^ M ^ mi + -f^ mx + M 2 <5 m2 J + 
1 / M 3 M 2 M 2 



A/f3 \ 

+-±6 3 ni +M 1 M 2 6 mi 6 m2 +M 3 6 3 ns +0(1/N 4 ) 
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3. MULTISCALE EXPANSION OF THE LPKdV EQUATION 

The lattice potential Korteweg-de Vries (lpKdV) equation is given by [12]: 

V = (p- q + u n , m+ i - u n+ x >m )(p + q- u n+1)Tn+1 + u ntm ) - (p 2 - q 2 ) = 0, (3.1) 

where p, q, p ^ q, are two real parameters. The above equation is probably the best-known 
completely discrete nonlinear equation which involves just four points which lay on two or- 
thogonal infinite lattices and it is nothing else but the nonlinear superposition formula for the 
Korteweg-de Vries equation. 

By defining fj, = p — q and ( = p + q, eq. (|3.ip can be written as 

V = [fi(T n T m u -u) + ({T n u - T m u)\ - [{T n u - T m u)(T n T m u - u)\ = Vi - V n i = 0, (3.2) 

where Vi and V n t denote respectively the linear and the nonlinear part of the lpKdV equation. 
The linear part Vi has a travelling wave solution of the form u = exp {i[nn — uj(K)m]} with 

uj(k) = —2 arctan ( - - tan — ) . (3-3) 

Applying the expansions (|2.7l 12.81 12.9|) to the function u = u ( m \ K > Pi takes the form 

OO . 

i=0 



where the operators L, = Li(T n ,T m , 5 ni , {5 mj }^ =1 ) can be constructed in a recursive way. The 



u «i j \ u m,j i 

lowest operators Li, i = 0,1,2, read 

Lq = n{T n T m -l) + ((T n -T m ), 

Li = ixT n T m ( M\8 ni + Mi5 mi ) + C ( MiT n 6 ni - M{I m b mY , . 



L% — fj,T n T m 



-^-8 2 ni + MiM^Sm, + -^-8 2 mi + M 2 5 m2 



+ 



+ c 



^Sl-^ ^$5 2 mi+ M 2 5 m2 



2 

As nonlinearity generate harmonics, let us now expand the function u as 

u = u {a) (m , {mi}f =l ; N y<-n~um) _ 

As u is assumed to be real then 

u (-a) = s (o) 

where by u we denote the complex conjugate of u. 
Moreover, if the nonlinear part should enter as a perturbation in the multiscale expansion of 
eq. (|3.2p . we need u( a )(rti, {mi}f =l ; oo) = 0, Va G Z. This implies that we have to expand each 
function vS°^ in inverse powers of N: 

oo 1 

« (Q) (ni,{m J }iI 1 ;iV) = ^-_4 a )(n 1 ,{m 4 }£ 1 ). 



N k 
fc=l 



Then reads 

i-l 



^ = EE 4 E L ^t\inx, {mi}f =l y^ n -" m \ 



N l 

aeZ i=l fc=0 



Performing the multiscale expansion of eq. (|3.2p we get several determining equations obtained 
selecting the different powers of 1/N and the different harmonics a. 

So, let us write down the resulting determining equations at the lower orders of 1/N for the 
harmonics a = 0, 1, 2, necessary to get a dNLS equation as a secularity condition. 
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The order 1/N gives, for a = 0, 1, linear equations which are identically satisfied by taking 
into account the dispersion relation (|3.3|) . For \a\ > 2 one gets some linear equations whose only 
solution is given by u[ a) = 0. 

The order 1/N 2 gives, for the harmonics a = 0, 1, 2, the following equations 

{H + C)Mi* ni + C)Mi<W] ^i 0) = 2 (-e iK - e" iK + + e^) | 2 , (3.4) 
e i« (^ e -i- + C ) M 1 5 ni v^ ) + e~' luJ {fie iK - () Mi6 mi v^ = 0, (3.5) 
[C (e 2iK - e~ 2iw ) + ii (e 2i ( K -^ - l)] 4* } = (-e iK + e"^ + e^"-") _ e i(*-2w)) ( n W)2_ (3 6 ) 

The solution of eq. (|3.5() is given by 4 ( n l> {m,j}^ 1 ) = 4 ( ra 2j {mi}f =2 ) with n2 = n\ + 7ml, 

7 = =pl, provided that the integers M\ and M\ are chosen as 

Mi = 7 5e" iw (/ie iK - C) , Mi = -5e iK (/ie"* 1 " + C) , (3.7) 

where 5 G C is a constant. As it has been shown in [9] one can always choose S = r exp (i#), 
with r > and 9 = — arctan [(£sin/-v)/(£cosK — //)], in such a way that M\ and Mi are indeed 
positive integers. Taking into account the dispersion relation (j3.3|) . the coefficients Mi and Mi 
in eq. (|3.7f) can be rewritten as 

M 1=1 S^-Ce iK ), Mi = 5e iK ^— 4. (3.8) 

Eqs. (|3.44I3.6| ) allow to express and vSp in terms of 4^ and 4 . As is a function of n 2 

the same must be for uf and i.e. «f (m, {m,}^) = 4 0) (n 2 , ^(m, W^) = 

u^\n2, {mi}f =2 ). Then, taking into account eq. (|3.8|) . we find that eqs. (|3.44I3.6|) reduce respec- 
tively to 

x ..(0)_ n , (i)| 2 „ ^ 2 7 (l + e iK ) 2 ,„ q x 

4 2) = a 2 rf } ) 2 , -^ d-^+o - (3 - 10) 

We can now consider the equation for the harmonic a = 1 at order 1/iV 3 . We have: 

(o"l^m + CT2^ mi ) + (cr 3 ^ a + cr 4 ^ i + fJ 5 (5 ni 5 mi + CT 6 5 m2 ) 4 = (3.11) 
i 1 ) r x 1 e x (0) . _(i) (2) 

where, taking into account eq. (13, 8D . the coefficients cj,, 1 < i < 9, read 



7 Se iK (// - Ce iK ) {p 2 - C 2 ) 

^ = J^— ( > ° 2 = -7*1. 

^^(M-C^K a4= 2(^-W) ' 
2 

°-B = - (M 2^2) » ct 6 = M2 (A* - Ce iK ) , 

. ai(e 2iK -l) . Se iK (/u 2 - C 2 ) (// + C) (1 - e 2iK ) 

a 7 = — 577 °8 ~ 



Ca* (e 2iK -l) 2 (e iK + l) 2 (Ai-C) 



<7 9 



e iK (/i - Ce iK ) (/ie iK - C) 2 
Using eqs. (13.91 13.X0H we can write eq. (13. lip as 



{o x 8 ni + cr 2 6 mi ) u 2 l) = C(u ( l y ), (3.12) 
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where C is a linear operator. Notice that the l.h.s. of eq. f)3. 12|) is the same as in formula (|3.5p . 
but it involves the field i4 > instead of . Requiring that no secular term appears we get the 
following equation for = u^\ri2, {Tni]f =2 ): 

(^ni + °"4^i + ^^i^mi + V&$m 2 ) «! = Uy ((7 7 5 ni + er 8 <W) «1 + <7g«i *4 • ( 3 - 13 ) 

Then from eq. f)3. 12|) we see that u 2 wu l satisfy the equation (<7i<5 ni +o" 2 <5 TOl )«2 = 0, i.e. 

1*2 ( n i>{ m i}^i) = ^2 ( n 2j {^i}^) whenever eq. f|3.8j) holds. Using eqs. (|3.9l 13.10]) we find 
that eq. f)3. 13j) is equivalent to the dNLS equation 

iS m2 uP = Pl S 2 n2 uP + P2 u^\uP\ 2 , (3.14) 

where 

i(pS 2 e iK (( 2 - /i 2 ) (e 2ire - 1) 



Pi 



P2 



2M 2 (^e iK - C) 
iCMM-C) (e 2iK -l) (e iK + l) 



M 2 e iK (/x + C) (ve iK - C) 2 (Ce iK - ' 
Taking into account the form of S, one finds that p\ and /? 2 are both real numbers: 

pO" 2 (C 2 - A* 2 ) sulk 



Pi 



P2 



M 2 (C 2 + P 2 - 2Cpcos «) 

8Cp(C — p)(l + c °s k) 2 sin k 



M 2 (p + C) (C 2 + ^ ~ 2CP cos k) 2 ' 

Notice that p\p 2 < 0, so that eq. (|3.14p is a defocusing dNLS equation. 

Eq. (|3.14|) may be obtained as a compatibility condition of the linear problem 

5 n2 <S>n2,m 2 (v) = U(uy ,uy;r))<f> m>m2 (r]), (3.15) 
Sm 2 ®n 2 ,m 2 (r]) = V (u^ , ^ ; r/)$ n2 , ma (ry) , (3.16) 



with 



/ ■ (i) 
,(i) „-(i).^ - I l7 l u i 

\ —irj 



. (i) _(i) . . / 2\7] 2 + \\ui ] \ 2 2r/«f ) - iS n2 u^ ) \ 

m«i»«i;»7) = „ _(i) ' _(i) 2 . f i) 2 , 

V 2r/u\ ' + i5 n2 n\ ; —2itj z - i\u\ ' \ z J 

being 77 G C the spectral parameter and <&n 2 ,m 2 (j]) a discrete complex vector function. A 
necessary and sufficient condition for the integrability of eq. (|3.14|) is that the coefficients p\ 
and p 2 be real. Thus eq. (|3.14|) is integrable iff i = 00. 

4. Concluding remarks 

In this paper we have extended the results obtained in [7, 9] to the case of functions of slow- 
varyness of order infinity To do so we have used the connection between shift operators and 
infinite series of differential operators (|2.2|) . Using these formulas we have been able to easily 
recover the formula (|1.1|) . usually proved by combinatorial techniques [5]. Moreover, in analogy 
with the continuous case we can introduce partial shift operators, which are expressed in terms 
of 5 operators, namely infinite series of difference operators. 

The multiscale expansion of the lpKdV equation in terms of harmonics implies, as a con- 
dition for the non existence of secular terms, a NLS equation written in terms of 5 operators. 
Choosing 5 n2 in eq. (|3,14p as a series of symmetric differences and fixing 1 = 2 the resulting 
dNLS equation is equivalent to the one presented in [7, 9] , see eq. (jl.2p . In this way we have 
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shown that for any finite order of slow-varyness the obtained dNLS equation will be local and 
no summation terms will ever appear. 

As shown in Section 3, eq. (|3. 14[) has a matrix Lax pair (|3.151 [37T6]) . which is also expressed 
in terms of 5 operators. Such Lax pair has a reduction to any finite order of slow-varyness, pro- 
vided that the associated wave functions $ n2i?Tl2 (7/) have the same finite order of slow-varyness. 
However in the case of finite £ the Lax pair (13.151 13.16P reduces to difference equations and their 
compatibility is no more satisfied by the slow- varying approximation of eq. (|3.14p . In fact the 
difference operators A do not satisfy the Leibniz rule. 

This work still leaves many open problems on which we are working at the moment. Let us 
just mention the analysis of the solutions of the lpKdV equation obtained from the exact ones of 
the reduced equation and their comparison with the ones obtained by carrying out a multiscale 
expansion of the continuous potential KdV equation; the reduction of symmetries of the lpKdV 
equation; the reduction of differential-difference equations both integrable and nonintegrable. 

It is still an open problem to understand the role played by an integrable discrete equation like 
the nonlocal discrete NLS equation introduced by Ablowitz and Ladik [1] from the isospectral 
compatibility of a discrete analogue of the Zakharov-Shabat spectral problem (13.151 l3~T6|) : 



1 



At 



2(Ax) 5 



u n +i, 



2u r 



n-1 

n 



A 



k,m I ~\~ 



(4.1) 



+ ( Un+l,m+l Y[ ^ fc > m ~~ ^ u n,m+l + ^n-l,m+l J 
\ i— — i / 



± - 



k=— oo 

Un,rn {Un,m'^ j n+l,m ~\~ Un,rn+l'^n+l,m-i-l^) ~i~ 



~i~ ^n,m+l {^n—l,m'^n,m ~\~ ^n- l,m+l^n,m+l ) ~\~ 



n-1 



fc=— oo 



+ 2|u rljm | U n -\.\^ m -\-\ | £ ^k,m 2 |l£n,m+l | U n —\^ m J^J A 
k=— oo 

n n— 1 



k,m 



k=—oo 



k=— oo 



with 



A 



k,m 



1 ± kfc,m+l| 



&k,m '^•k,m'^'k—\,m ^ , fc+l,m'^fc,m' 



1 ± \uk,m\ 2 

Eq. ()4.ip reduces to the continuous NLS equation when At — > and Ax — > 
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